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Abstract 

A field-theoretic description of the critical behaviour of systems with quenched defects 
obeying a power law correlations ~ |x|~ a for large separations x is given. Directly for 
three-dimensional systems and different values of correlation parameter 2 < a < 3 a 
renormalization analysis of scaling function in the two- loop approximation is carried out, 
and the fixed points corresponding to stability of the various types of critical behaviour 
are identified. The obtained results essentially differ from results evaluated by double e, 5 
- expansion. The critical exponents in the two-loop approximation are calculated with 
the use of the Pade-Borel summation technique. 



PACS: 64.60.Ak; 64.60.-i; 64.60.Fr 



In last years, much effort has been devoted to investigation of the critical behaviour of 
solids containing quenched defects. In most papers considerations have been restricted to 
the case of point defects with small concentrations so that the defects and corresponding 
random fields have been assumed to be Gaussian distributed and 5 - correlated. 

For the first time in work of Weinrib and Halperin (WH) |Ij have been offered the 
model of the critical behaviour of disordered system in which the correlation function of 
the random local transition temperature g(x — y) =<< T c (x)T c (y) >> — << T c (x) >> 2 
falls off with distance as a power law ~ |x — y| _a . It was shown that for a > d long-range 
correlations are irrelevant and the usual short-range Harris criterion 2 — dv Q = a Q > 
of influence of 5 - correlated point defects is realized, where d is the spatial dimension, v 
and a Q are the correlation-length and the specific-heat exponents of the pure system. For 
a < d it was established the extended criterion 2 — av Q > of disorder influence on the 
critical behaviour. As a result, a wider class of disordered systems, but not only three- 
dimensional Ising model with 5-correlated point defects, can be characterized by new type 
of critical behaviour. So, for a < d it was discovered a new long-range (LR) disorder stable 
fixed point (FP) of the renormalization group recursion relations for systems with number 
of components of the order parameter m > 2. The critical exponents were calculated in 
the one-loop approximation using a double expansion in e = 4 — d <^il and 5 = 4 — a <C 1 . 
In the case m = 1 the accidental degeneracy of the recursion relations in the one-loop 
approximation did not permit to find LR disorder stable FP, but it was predicted for 
5 > 5 C = 2 (6^/53) ly/2 a change in critical behaviour of the model from short-range (SR) 
to LR correlation type. Korzhenevskii, Luzhkov and Schirmacher |§ have proved the 
existence of the LR disorder stable FP for the one-component WH model and found the 
characteristics of this type of critical behaviour. Also they have considered very interesting 
model of the critical behaviour of crystals with LR correlations caused by point defects 
with degenerate internal degrees of freedom ||, f|. 

The models with LR correlated quenched defects present both doubtless theoretical 
interest from possibility of prediction a new types of the critical behaviour in disordered 
systems and experimental interest from possibility of realization RL correlated defects in 
the orientational glasses || and disordered solids containing defects of fractal-like type 
However, numerous investigations of pure and disordered systems performed with 
the use of the field-theoretic approach show that the predictions made in the one-loop 
approximation, especially on the basis of the e - expansion, can differ strongly from the 
real critical behaviour [6-9]. Therefore, the map of regions with the various types of 
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critical behaviour received for WH model on the basis of e, 5 - expansion [|TJ (figure 1(a)) 
may be not corresponding to the critical behaviour of the three-dimensional WH model 
for different values of m and a. In this case the results for the models with LR correlated 
defects received with use of e, 8 - expansion [1,3,4,10-12] must be corrected. To shed 
light on this question and to determine more accurately the dependence of the critical 
behaviour on the number of components of the order parameter m and the values of 
correlation parameter a, we have constructed a field-theoretical description of the three- 
dimensional WH model in the two-loop approximation for the values of a in the range 
2 < a < 3. 

The effective Hamiltonian of WH model after using the replica trick is given by 
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where 4> a = YhL\ §i a -, 4>ia is (n x m)-component order parameter. The properties of the 
original disordered system are obtained in the replica number limit n — > 0. The Fourier 
transformation of the interaction vertex g(x) ~ x~ a gives g(k) = v + wok a ~ d for small k. 
g(k) must be positive definite, therefore if a > d, then the w term is irrelevant, v o > 
and H ef[ (1) corresponds to model with SR - correlated defects, while if a < d, then the 
w term is dominant at small k and wq > 0. 

As is known, in the field-theoretic approach [|1^] the asymptotic critical behaviour of 
systems in the fluctuation region are determined by the Callan-Symanzik renormalization- 
group equation for the vertex parts of the irreducible Green's functions. To calculate the 
(3 functions and the critical exponents as functions of the renormalized interaction vertices 
u, v and w (scaling 7 functions) appearing in the renormalization-group equation, we used 
the standard method based on the Feynmann diagram technique and the renormalization 



procedure ||14j| . The three types of interactions can be represented graphically as in figure 
2(a). When we considered a diagrammatic representation of two-point vertex function 
V( 2 \ three types of four-point vertex functions and two-point with the 2 insertion 
vertex function r^ 1 ' 2 ) in the two-loop approximation the diagrams were integrated numer- 
ically in d — 3 and with values of parameter a determining momentum dependence of the 
w interaction in the range 2 < a < 3 with changes through the step Aa = 0.01. Unlike 
the works using e, 5 - expansion we took into consideration the graphs of the form (figure 
2(b)), contributions of which are increased when the values a are removed from a = 3. 

As a result, we obtained the (3 and 7 functions in the two-loop approximation in the 
form of the expansion series in renormalized vertices u, v and w. Because of impossibility 
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in short notes to present the coefficients of these series for different values of a we give 
here the obtained f3 and 7 functions only for a = 2 (the case with a = 2 corresponds to 
system of straight lines of impurities or straight dislocation lines of random orientation 
in a sample): 



f3 u (u, v, w) = -u + u 2 - \uv - 1.901416m, - fggpV + ^fgf 1 ^- 

_ ig uv2 + (1.230378^+6.713002) A _ 031265W _ 1.19 347W) 

P v (u, v,w)=v + v 2 + § w 2 + 1.901416™ - + ||f 3 + 0.488229w 3 - 

-f^r§^ 2 - 1.974883^±|W + ^4m 2 ^ + 0.806375^ 2 + 

27(m+8) (m+8) 27(m+8)' ! 



+0.839125t» 2 w 



; (m+2) 
(m+8) ' 

A>, v, w) = 2w; + 0.628176W 2 + \vw - 2 -^uw - 0.1528u; 3 + §^±^u 2 w+ (2) 
+ ^^ 2 w + 0.090516ggww 2 - 0.022629™ 2 - §^uvw, 

l4> (u, v, w) = 0.004222w + ^ + + % u 2 + ^v 2 + 0.056893w 2 - $0^ U v- 
-0.315823^±§ww + 0.078956m/;, 

(m+8) ' 

7 ^(u, v, w) = -%&u +\v + 0.31831^ + j£0u 2 + ±v 2 - 0.019507^ 2 - 

~ °- 270565 TSIy^ + 0.067641™. 

The series (2) are normalized by a standard change of variables pi M u—>-Qu/('m + 8)J, 
v — > v/32J, w — > w/32J, so that the coefficients of the terms u, u 2 and v, v 2 in f3 u and /3„ 
become 1 in modulus, where J = J d d q/(q 2 + l) 2 is the one-loop integral. 

The nature of the critical behaviour is determined by the existence of a stable FP 
satisfying the system of equations 

Pi(u*,v*,w*) = (1 = 1,2,3). (3) 

It is well known, that perturbation series are asymptotically convergent, and the vertices 
describing the interaction of the order parameter fluctuations in the fluctuating region 
r — > are large enough so that expressions (2) can't be used directly. For this reason, to 
extract the required physical information from the obtained expressions, we employed the 
Pade-Borel approximation of summation of asymptotically convergent series extended 



to the multiparameter case ||, [15]] . We used the [2/1] approximant to calculate the (3 
functions in the two-loop approximation. 

However, the analysis of the series coefficients for (3 W function has shown that the 
summation of this series is fairly poor, which resulted in absence of FP with w* 7^ 0, for 
example, in the case m = 1 for a < 2.93, in the case m = 2 for a < 2.67 and so on. 
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Dorogovtsev has found the symmetry of the scaling function for WH model in relation to 
the transformation (u, v,w) — > (u,v,v + w) [K| which gives the possibility to investigate 
the problem of FP existence with w* 7^ in the variables (u,v,v + w). In this case our 
investigations have shown the existence of FP with w* 7^ in the whole region where the 
parameter a changes. 

We have found three types of FP's in the physical region of parameter space u*, v*, v* + 
w* > for different values of m and a. The type I corresponds to FP of pure system 
(u* ^ 0,v*,w* = 0), the type II is SR-disorder FP {u\v* ^ 0,w* = 0) and the type III 
corresponds to LR-disorder FP's (it*, v*, w* 7^ 0). The type of the critical behaviour of this 
disordered system for each value of m and a is determined by stability of corresponding 
FP. The requirement that the fixed point be stable reduces to the condition that the 
eigenvalues of the matrix 

dPi{u\,u* 2 ,ul) 

= — ^ — (4) 

lie in the right-hand complex half plane. 

Values of the stable FP's obtained for the most interesting values of the number of 
order-parameter components m and 2 < a < 3 are presented in table 1. As one can see 
from this table, for Ising model (m = 1) the LR-disorder FP is stable for values of a in 
the whole investigated range. The additional calculations for 3 < a < 4 have shown that 
only FP II is stable in this range. For a = 3 FP values for vertices u and g(k) are equal 
u* = 2.383 38, g* = v* + w* = 0.55164 and correspond to SR-disordered Ising model 
FP, although w* 7^ 0. Like that, for m = 1 and a = 3 the LR disorder is marginal, and 
the critical behaviour of WH model, as of SR-disordered Ising model, is characterized by 
the same critical exponents (table 2). The critical behaviour of the XY- model (m = 2) is 
determined by the LR-disorder FP for a < 2.96 and the SR-disorder FP for a > 2.96. The 
Heisenberg model (m = 3) is characterized by the change of the critical behaviour types 
from the LR-disorder type (III) for a < 2.85 to the pure type (I) for a > 2.85. Figure 
1(b) shows regions of the various types of critical behaviour of the WH model, which we 
obtained in the two-loop approximation. The large change in the picture indicates that the 
correspondence between WH results and our calculations in the two-loop approximation 
is weak. 

However, the results, which we received for disordered XY-model, must be corrected. 
We think that in the higher field-theory orders of approximation k the critical behaviour 
of XY-model will be determined by the FP of pure type (I) for a[ k ^ < a, but not by the 
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SR-disorder FP (II), obtaining in the two-loop order. Here, a^ is marginal value for a 
in k-th. order of approximation, for which disorder is irrelevant (a^ ~ 2/v Q = 2.99 with 
v a = 0.669 |L(J for m = 2). The two facts indicate this, such as the weak stability of 
the SR-disorder FP revealed for 2.96 < a < 4 and that the = 3 for m c = 2.0114. In 
the higher orders of approximation the marginal value of m c can be found with the use 



of the Harris criterion dv (m c ) — 2 = 0, and such as u = 0.669 [|1J for m — 2, then 
m c < 2. Therefore, we think that the corrected picture of the regions of various types 
of the critical behaviour of the model with LR-correlated defects will be represented by 
figure 1(c). 

Finally, we have calculated the static critical exponents for the WH model (table 2), 
received from the resummed by the generalized Pade-Borel method 7 functions in the 
corresponding stable FP's r] = j^u* , v* , w*) , v = [2 + 7^2(1**, v*, w*) - ^(u*, v*, w*)]^ 1 . 
Also, we have found by the method of the work M the dynamic scaling function 7a and 
calculated the values of the dynamic exponent z = 2 + 7a(w*, v *, w*) on the basis of the 
resummed 7a function (table 2). As example, for a = 2 the received 7a function is given 
by 

7a(u, v, w) = \v + 0.314088w + 0.226777g^w 2 + §^v 2 - 0.0764w 2 - 
-0.092593 { i^uv + 0.123604^±|ww - 0.011315mo. 

(m+8) (m+8) 

The comparison of the exponent v values and ratio 2/a from table 2 shows the violation 
of supposed in |I| on the basis of some heuristic arguments as exact the relation v — 2/a. 
The revealed difference is caused by the use in our work of more accurate field-theoretic 
description in the higher order of approximation for three-dimensional system directly 
together with methods of series summation. Also, these distinctions can be explained by 
the application for calculations of the concrete numerical values of parameter a and the 
taking into consideration the graphs of the form (figure 2(b)), thrown away when e, 5 - 
expansion is used, but contributions of which are increased when the values a are removed 
from a = 3. Of course, there are errors of the present consideration determined by the 
accuracy of series summation for the (3 and 7 functions. However, the comparison of 
exponents values for SR-disorder Ising model, calculated with the use Pade-Borel method 
in 0] and in the two-loop and four-loop approximations accordingly, shows that their 
differences not more than 0.02. In the same time, in our work v — 2/a depends on the 
values of a and m and has the value 0.284, as example, for a = 2 and m — 1, that is 
considerably larger. 

In closing, we hope that the features of the critical behaviour of the WH model revealed 
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in our paper will stimulate the organization of experimental works in real disordered 
systems with long-range correlated defects like the orientational glasses and solids with 
defects of fractal-like type. Also, the computational methods can be applied for simulating 
disordered systems with straight lines of impurities of random orientation in a sample 
(a = 2). The received values of exponents can be used for explanation of the results 
of computer simulation of the three-dimensional disordered Ising model |I7j at impurity 
concentrations between the threshold of impurity percolation and the spin-percolation 
threshold, in which the fractal-like behaviour of impurity extended structures and the 
competition between impurity-percolating and spin-percolating clasters are possible. 
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Figure and table captions 



Figure 1: 

Regions of the various types of critical behaviour, which have been determined (a) in 
U] on the basis of the double e, S - expansion; (b) in the present paper with use of the field 
theoretic description in two-loop approximation for three-dimensional WH model; (c) in 
the present paper with taking into consideration of the higher orders of approximation. 

Figure 2: 

Graphs () that correspond to vertices u, v and w ; (b) that are taken into consideration 
in addition as compared with works, using e, 5 - expansion, H corresponds to vertices 
u, v and w. 

Table 1: 

Stable fixed points of the 3D WH model from two-loop expansions. 
Table 2: 

Critical exponents of the 3D WH model from two-loop expansions. 
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Table 2 
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